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The B° s — > J/ipK + K channel has previously been studied only when the K + K are 
consistent with the decay of the meson. This mode has been used to measure the CP 
violation in B s mixing, S , a key probe in the search for physics beyond the Standard 
Model [H 121 |3l H] Q In addition to the other resonant or non-resonant final states may 
appear and affect the CP measurements, including an S-wave contribution |5j. In this 
paper we study the entire K + K~ mass spectrum, including a search for other final states 
that may be useful in CP violation studies. These states may provide other ways of 
measuring <j) s , in decays with a different spin structure that may be useful for revealing 
different aspects of CP violation. 

We use a O.lGftr 1 data sample collected with the LHCb detector [UJ at a center-of- 
mass energy of 7 TeV. The detector elements most important for this analysis include 
a vertex locator, a silicon strip device that surrounds the pp interaction region in the 
LHC, and other downstream tracking devices before and after a 4Tm dipole magnet. 
Two ring-imaging Cherenkov detectors are used to identify charged hadrons, while muons 
are identified using their penetration through iron. This analysis is restricted to events 
accepted by a di-muon trigger [7]. Subsequent selection criteria are applied that serve 
to reject background, yet preserve high efficiencies, as determined by Monte Carlo (MC) 
events generated using Pythia [8J, and the LHCb detector simulation based on Geant 
[9]. To be considered as a J/ip — > candidate, opposite sign tracks are required to 

have transverse momentum, p^, greater than 500 MeV, be identified as muons, and give 
a good fit to a common vertex]^] Di-muon candidates with masses between —48 and +43 
MeV of the J/ip peak are selected for further analysis, where the r.m.s. resolution is 13.4 
MeV. The invariant mass of the pair is constrained to the J/ip mass for further 

analysis. 

Kaon candidates are selected if their minimum distance from the closest primary vertex 
is inconsistent with being produced at that vertex. They must be positively identified 
based on the logarithm of the likelihood ratio comparing two particle hypotheses (DLL). 
There are two criteria used; loose corresponds to DLL(i^ — 7r) > 0, while tight has 
DLL(fC — 7r) > 10 and DLL(K — p) > —3. We use the loose criterion for checking kaon 
identification efficiencies, otherwise the tight criterion is used. In addition, the two kaons 
must have the sum of the magnitudes of their p^ > 900 MeV. 

To select B s candidates we require that the K + K~ pair and the J/ip candidate give 
a good fit to a common secondary vertex with a x 2 < 5 per degree of freedom. We also 
require that the B s candidate's decay point must be > 1.5 mm from the primary vertex 
and that the negative of its momentum vector points back to the primary. 

The B s candidate invariant mass is shown in Fig. 1 A clear signal is seen, part of 
which comes from the previously known J/ip<p> mode. A check was made for possible 
resonant states decaying to J/ipK~ since similar exotic states have been claimed [TU] . 
but no obvious structures are visible in the invariant mass spectrum. Figure [2] shows 
the K + K~ invariant mass for both signal and sideband regions, where the signal region 



1 Charge conjugate modes are also considered throughout. 
2 We work in units where c = 1. 
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Figure 1: Invariant mass of J/ipK + K combinations. The vertical lines indicate the 
signal and sideband regions. 



extends ±25 MeV around the B s mass peak and the sidebands extend from 35 MeV to 60 
MeV on either side of the peak. Apart from the large peak at the <fi there is a structure 
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Figure 2: Invariant mass of K + K combinations. The histogram shows the data in the 
signal region while the points (red) show the sidebands. 



near 1525 MeV. In addition there is an excess of signal events over most of the mass range 
which we refer to as non-resonant. We investigate the possibility of the peak to be the 
/2(1525) resonance. The PDG quotes the mass of the f 2 state as 1525±5 MeV and the 
width as 73±g MeV [TT| - Other states such as the / 2 (1270) and the / (1500) have small 
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branching fractions into K + K of less than 5%, and are unlikely to have large rates. 

It is possible for the decay B — > J/iIjK~7t + to fake our signal if the 7r + is misidentified 

as a K + . A specific example is given by B — > J /ipK 2 (lA'$Q) decays [12]. To examine 

if we are sensitive to a reflection of this mode in the 1525 MeV di-kaon mass region, 

a simulation was performed where the 7r + from the K 2 (14:30) was interpreted as a K + . 

Figure [31(a) shows that the reflection of this mode does indeed peak in the di-kaon mass 

— o 

range around 1525 MeV. It also peaks in the B s signal region but is much wider than the 
B s mass peak. The region 25 — 200 MeV above the B s mass peak provides a sample of 
misidentified B — > J/tfjK~ir + decays. By measuring the number of B candidates in this 
region we can calculate the number in the B s signal region. 




m(K + K') (MeV) m(J/\|/K"re + ) (MeV) 



Figure 3: (a) The m(K + K ) distribution for simulated B — > J/ipK* 2 ( 1430) decays where 
the 7r + from the i^ 2 (1430) decay is interpreted as a K + . The histogram shows m(K + K~) 
in the signal region of B s mass and the points in the sideband region. The simulation 
corresponds to approximately 8 fb _1 of data, (b) The m(J/if)K + 7c~) distribution for 

J/ipK + K~ data candidates 25 — 200 MeV above the B s mass, and with m(K + K~) within 
±300 MeV of 1525 MeV, reinterpreted as B -»■ J/ipK-n + events. The fit is to a signal 
Gaussian whose mass and width are allowed to vary as well as a quadratic background. 



To determine the size of any B° reflection in the f 2 mass region we select events where 
the reconstructed J/ipK + K~ mass is in the range 25 — 200 MeV above the B s mass, 



reassign each of the two kaons in turn to the pion hypothesis, and plot the J '/ipKir mass. 

The resulting peak at the B mass has 36±10 events from the fit shown in Fig. 3[b). 

— o 

We calculate 37±10 events in the B s signal region, using the shape from Monte Carlo 
simulation, and use this number as a constraint in the fit described below to determine 
the f 2 parameters and signal yields. 

To test the f 2 hypothesis we perform a simultaneous fit to the B s candidate mass and 
the di-kaon mass. The f 2 signal is parameterized by a spin-2 Breit-Wigner function [T5] . 
The width of the f 2 is fixed to the PDG value of 73 MeV A contribution from non- 
resonant K + K~ is included as a linear function in the di-kaon mass. The contribution 
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Figure 4: Projections of fits to (a) the B s candidate mass and (b) the di-kaon mass. The 
f 2 signal is parameterized by a spin-2 Breit-Wigner function whose width is fixed to 73 
MeV (dotted curve). The combinatorial background is shown in the light shaded region, 
while the darker shaded region shows the non-resonant J/ipK + K~ component. The long- 
dashed (red) line shows the misidentified B — > J/iI)K~-k + decays, and the (blue) line the 
total. 



from the K~ti + reflection is included using the di-kaon and B s mass shapes from the 
simulation, with the normalization fixed by the event yield in Fig[3|b). The results of the 
fits are shown in Fig. |4j The f 2 mass from the fit is 1525±4 MeV and the yield is 269±26 
events within ±125 MeV of the mass. If we allow the f 2 width to vary we find a consistent 
value of 90^4 MeV. As we have not taken into account possible interferences between the 
f 2 and other J/ipK + K~ final states we do not provide systematic uncertainties for these 
values. The decay angle of the J/ip, Oj/^, can test for pure spin-0, or the presence of a 
higher spin state such as the spin-2 f 2 [IT] . Here 9j/^ is defined as the angle of the /i + 

with respect to the B° s direction in the J/ip rest frame. It is distributed as 

/(cos 6 J/4,) = {l-p) sin 2 6j/tp + - (l + cos 2 9j/^j , (1) 

where 1 — p is the fraction of helicity zero and p is the fraction of helicity ±1. Shown in 
Fig. [5] is the background subtracted, acceptance corrected cosOj/^, distribution for K + K~ 
masses in the f 2 region. MC simulation is used to find the acceptance correction. The 
points are extracted from the joint fit to the m(J/if)K + K~) and m(K + K~) distributions 
in the K + K~ mass region within 1400 — 1650 MeV for events in the peak above the 
non-resonant K + K~ . The fit result is p = (0.57 ± 0.13), with x 2 /number of degrees of 
freedom (ndof) of 10/8 (27% probability). Fitting only with an S-wave gives xV n dof of 
27/9 (0.1% probability), showing that the data are not likely to be pure spin-0, but are 
compatible with a higher spin state consistent with an f 2 contribution. 

The branching fraction of B s — > J/ipf 2 relative to B s — > J/ip<p is determined by 
assuming that the dominant background is S-wave and the signal D-wave, so there is no 
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Figure 5: Distribution of cos 9 j/^p for B° s — > J/iftfy decays. The background and non- 
resonant K + K~ components have been subtracted, and the data have been corrected for 
acceptance. The fit to Eq. [T] is shown by the solid line. Note that for pure S-wave the 
distribution would be sin 2 9 (p — 0), shown as the dotted curve, while for pure helicity 
1 (p — 1) the data would be described by the dot-dashed curve. 



interference between themfl The number of J/i/jK + K events is determined by a fit to 
— o 

the B s mass distribution, within ±20 MeV of the <p mass. A small S-wave component 
in the (f) mass region of (4.2±2.3)% is subtracted [2]. Although there are the same final 
state particles in both modes, the relative efficiency is (78±2)%, where the uncertainty 
arises from simulation statistics. The efficiency ratio differs from unity due to the different 
Pt distributions of the kaons in the final states. The kaon identification efficiencies are 
corrected with respect to those given by the MC simulation using a sample of D* + decays, 
where the kaon can be selected without resorting to PID information. Typical corrections 
are on the order of 5%. 

To find the effective relative rate of f 2 decays we use the fit where the width is allowed 
to vary. There are 320±33 f 2 events and 1774±42 <fi events. Correcting for the relative 
efficiencies and the explicit branching fractions B (f 2(1525) — > K + K~) = (AAA ± 1.1)%, 
and B (<p ->• K+K~) = (48.9 ± 0.5)% [EE], we measure 

B (b° s J/^(1525)) 

R = V r — '- = (26.4 ±2.7 ±2.4)%. (2) 

B (W s -> 

The systematic uncertainty on R has several contributions, as listed in Table [1} The 
largest source of uncertainty is f 2 width. The error quoted reflects changing the width 

3 Although there can be interference as a function of the K + decay angle in the rest frame, integrating 
over this variable causes the net result to be zero. 
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by one standard deviation from the fitted value of 90 MeV. The helicity amplitudes of 
the J/ipf'2 decay are unknown, unlike the J/i[>4> amplitudes which are well measured [TTj . 
The difference between the values obtained using helicity zero and helicity one J/ip MC 
samples is 4% compared to our central value. The S-wave subtraction of the events in the 
J/ijxp region causes a 2.3% uncertainty. We include an uncertainty for the efficiency as 
a function of K + K~ mass, as the tracking could be sensitive to the opening angle of the 
kaon pair. Modifying the acceptance from a flat to linear function of mass changes the 
yield by 2.3%. Varying the B s distribution within limits imposed by the data results 
in a small 0.5% change in the rate. Changing the mass resolution by its error results in 
a 0.5% change. A PID uncertainty of 1% is added to account for different momentum 
distributions of the kaons in the two final states. As a check we note that the ratio of 
the number of events in J/ifxp with tight cuts to loose cuts on the kaon identification is 
(61±2)% and the simulation gives a consistent (60±1)%. Variation of the background 
and signal shapes makes small differences. 

In conclusion, we have made the first investigation of the B° s — > J/if)K + K~ final state 
over the entire range of K + K~ mass. There is a significant non-resonant component that 
extends under the region which can affect CP violation measurements [5] . We have also 
observed B s — > J / ip ^(1525) decays. The branching fraction ratio relative to J/ip4> is 

B (b° s J/^(1525)) 

V ,_ n = (26.4 ±2.7 ±2.4)%, (3) 

B [B° s -> J/^J 

assuming that the background does not interfere with the signal amplitude. This decay 
mode can also be used to measure CP violation in the B s system, although a different 



Table 1: Systematic uncertainties on R. 



Source 


Change (%) 


f 2 width 


6.3 


Helicity 


4.0 


Relative efficiency 


2.6 


S-wave under <p 


2.3 


K + K~ mass dependent efficiency 


2.3 


Background shape 


1.3 


B s p T distribution 


0.5 


B s mass resolution 


0.5 


PID 


1.0 


Signal shape 


1.0 




2.5 


B{<j>^ K + K~) 


1.0 


Total 


9.2 
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transversity analysis than in J/ijjcj) would be required as the final state is a combination 
of a spin-1 J/ip and a spin-2 f 2 state. Some consideration has been given to measuring 
CP violation in vector-tensor decays [H]. 
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